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Active fluids such as swarming bacteria and motile colloids exhibit exotic properties different 
from conventional equilibrium materials. As a peculiar example, a spherical tracer immersed inside 
active fluids shows an enhanced translational diffusion, orders of magnitude stronger than its intrin¬ 
sic Brownian motion. Here, rather than spherical tracers, we investigate the diffusion of isolated 
ellipsoids in a quasi-two-dimensional bacterial bath. Our study shows a nonlinear enhancement of 
both translational and rotational diffusions of ellipsoids. More importantly, we uncover an anoma¬ 
lous coupling between particles’ translation and rotation that is strictly prohibited in Brownian 
diffusion. The coupling reveals a counterintuitive anisotropic particle diffusion, where an ellipsoid 
diffuses fastest along its minor axis in its body frame. Combining experiments with theoretical 
modeling, we show that such an anomalous diffusive behavior arises from the generic straining flow 
of swimming bacteria. Our work illustrates an unexpected feature of active fluids and deepens our 
understanding of transport processes in microbiological systems. 

PACS numbers: 87.17Jj, 05.40-a, 47.63.Gd 


The diffusion of a small tracer in a surrounding 
medium provides a reliable means for probing material 
properties of complex fluids [1]. Particularly, such a 
method has been used to investigate unique features of 
active fluids—a novel class of nonequilibrium soft ma¬ 
terials with examples across a wide range of biological 
and physical systems including flocking animals |2-4[, 
vibrated granular beds [B|, synthetic colloidal swimmers 
H, 0, and a self-propelled cytoskeleton [8, 9]. The behav¬ 
ior of spherical tracers in active fluids is most clearly illus¬ 
trated by the diffusion of colloidal spheres in suspensions 
of swimming microorganisms 10l-ll8| . Because of the hy¬ 
drodynamic and steric interactions ghh , a colloidal 
particle immersed in a bath of microswimmers exhibits 
a super-diffusive behavior at short times and a dramati¬ 
cally enhanced translational diffusion at long times. 

The enhanced diffusion of passive particles such as 
nutrient granules, dead bacterial bodies and extracel¬ 
lular products is of great biological importance, which 
maintains an active ecological balance 0 , stimulates 
bio-mixing 0 , and promotes intercellular signaling and 
metabolite transports [23[. However, few natural par¬ 
ticles have the perfect spherical symmetry and usually 
possess more than translational degrees of freedom. It 
is still an open question how and to what extent the en¬ 
hanced translational diffusion of an anisotropic particle 
is influenced by other degrees of freedom especially by 
its rotation. Furthermore, as one of the most prominent 
features of active fluids, the emergent collective motion 
of self-driven units leads to swarming patterns 00, 
which strongly affect systems’ rotational degrees of free¬ 
dom. Thus, the rotational dynamics of anisotropic par¬ 
ticles also provide a way for probing intrinsic properties 
of the swarming phase of active fluids [26j]. 

In our experiments, we use suspensions of Escherichia 



FIG. 1: (Color online). Diffusion of an ellipsoid in quasi- 
two-dimensional bacterial bath, (a) Velocity field of swarm¬ 
ing bacteria around an ellipsoid obtained from particle image 
velocimetry (PIV). Scale bar: 20 /an. (b) Trajectory of an 
ellipsoid in a time interval At = 0.2 s. The trajectory is 
obtained from a custom algorithm for particle tracking ve¬ 
locimetry (PTV). The center of mass of the particle at differ¬ 
ent times, x{t), is indicated by black dots. The color indicates 
the orientation of the particle, 0{t), with respect to the x axis 
fixed in the lab frame. Bacterial concentration n = 30no- 


coli as our model active fluids (see supplementary mate¬ 
rial 0)- A bacterial suspension of concentration n is 
suspended in a free-standing soap film, where n is mea¬ 
sured in the unit of no = 8 x 10 8 cells/ml. The film is 
adjusted to be 15 =t 5 /am in thickness and 5x5 mm 2 in 
area. Polystyrene ellipsoids are used as our anisotropic 
tracer particles. In our study, we mainly focus on ellip¬ 
soids with semiprincipal axes of length a = b = 2.8 ± 0.2 
/am and c = 14.2T0.5 /am and aspect ratio p = c/a = 5.1, 
although two other ellipsoids with different sizes and as¬ 
pect ratios are also tested (see Ref. [I?} and Fig. EJ). 

Since the major axis of ellipsoids is larger than 
the thickness of film, particles undertake a quasi-two- 
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FIG. 2: (Color online). Enhanced translational and rotational diffusions, (a) Translational mean-squared displacements 
(MSDt) and (b) Rotational mean-squared displacements (MSDr). The slopes of lines are indicated, (c) Translational and 
rotational diffusion coefficients, Dt and Dr , as a function of bacterial concentrations, n. Solid lines indicate linear and nonlinear 
enhancements. The vertical line indicates the onset of obvious bacterial swarming. 


dimensional motion in the bacterial bath. The diffu¬ 
sion of ellipsoids is recorded with high-speed optical mi¬ 
croscopy (Fig.Hk), which allows us to extract the center- 
of-mass position, x(t) = [x(t),y(t)\, and the orientation, 
#(£), of ellipsoids at different times, t (Fig. Ud). 


First, we analyze the translation and rotation of el¬ 
lipsoids in the lab frame. Figures [2^,b show ellipsoids’ 
translational mean-squared displacements (MSDt), 
([A x(t)]) 2 = ([x(t + to) — x(to)] 2 ), and rotational mean- 
squared displacements (MSDr), ([A 0(t)] 2 ) = ([0(t+to) — 
6(to)] 2 } : at different n, respectively. The two degrees of 
freedom exhibit qualitatively similar trends. At short 
times, the particle motion appears to be super-diffusive 
following ([A x(t)] 2 ) ~ t ST and ([A 0(t)] 2 ) ~ t Sn with 
the power exponents 1.5 < Sr < 2.0 and 1.3 < Sr < 
1.6. At long times, the motion becomes diffusive with 
St,r ~ 1.0. The corresponding diffusion coefficients, de¬ 
fined in the diffusive regime as Dr = (\Ax(t)] 2 ) / At and 
Dr = ([A0(t)] 2 )/2t: increase linearly with n at low n 
(Fig.Efc), similar to the translational diffusion of spherical 
particles in dilute suspensions of microswimmers HM1. 
ffowever, for n > 10no, where bacteria develop collective 
swarming motions (Fig. [lji), the enhancement of Dr,r 
becomes nonlinear following a trend consistent with a 
power-law scaling Dr,r ~ n 2 (Fig. Et). 

The qualitative change of Dt,r(ti) scaling reflects the 
underlying transition of bacterial susp ensions from the 
disordered to the swarming phase [2J, |25| . Generally, a 
diffusion coefficient, Dr,r ~ / 2 , where l is the step size 
of particle motion in a short time in the superdiffusive 
regime. I is determined by the strength of the unbal¬ 
anced fluid flow created by numerous swimming bacteria 
at the location of ellipsoids, which in turn is proportional 
to the fluctuation of bacterial concentration surrounding 
the ellipsoids, Sn. Thus, Dr,r ~ l 2 ~ Sn 2 . The linear 
scaling Dr,r ~ n at low n indicates Sn ~ n , consis¬ 
tent with the central limit theorem In contrast, the 
nonlinear enhancement Drr ~ n 2 at high n indicates a 


giant number fluctuation Sn ~ n, which is known as a 
hallmark of the swarming phase of active fluids Slffl. 

Next, we investigate the coupling between the trans¬ 
lational and rotational motions of ellipsoids by measur¬ 
ing a dimensionless cross-correlation function, C(t) = 
(AxAy sin 20)/([Ax(t)] 2 } (Fig. [3^). C(t) deviates from 
zero, revealing a coupling between the two degrees of free¬ 
dom. More interestingly, C(t) shows a transition from 
positive to negative correlations with increasing n. In 
sharp contrast, the translation and rotation of a Brow¬ 
nian ellipsoid are always positively coupled [29[. The 
translational friction coefficient parallel to the major axis 
of an ellipsoid, £||, is invariably smaller than that per¬ 
pendicular to its major axis, (j_—a generic feature of 
Stokes flow [30]. Thus, from the Einstein relation, a 
Brownian ellipsoid diffuses faster along its major axis 
with D\\/D±_ > 1, where Di = ksT/Q is the diffusion 
coefficient in the body frame of the particle with i = ||j 
(parallel to the major axis) or i_ (perpendicular to the 
major axis) and ksT is the thermal energy. Hence, an 
ellipsoid undertaking Brownian motion has a tendency 
to diffuse along its major axis, resulting in a positive 
translation-rotation coupling in the lab frame as shown 
schematically by the blue ellipsoid in the inset of Fig. [3 }d. 
Note that when diffusing faster along its major axis, the 
ellipsoid rotates counterclockwise, leading to a positive 
0 and, therefore, (Ax Ay sin 2$) > 0. In contrast, the 
negative coupling indicates D\\/D±_ < 1 (the red ellip¬ 
soid in Fig. [3 b inset), a counterintuitive result strictly 
prohibited in Brownian diffusion. Thus, anisotropic par¬ 
ticles in bacterial suspensions possess an unusual mode 
of translation-rotation coupling, non-existent in equilib¬ 
rium systems. 

Motivated by the above analysis, we directly mea¬ 
sure anisotropic particle diffusions in the particles’ body 
frame. We transform the trajectories of ellipsoids into the 
body frame and obtain MSDs parallel and perpendicular 
to their major axes [13 • The body-frame MSDs show the 
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FIG. 3: (Color online). Translation-rotation coupling, (a) Cross correlations at different n. The dashed line indicates zero 
coupling, (b) The ratio of diffusion coefficients along the major and minor axes in the body frame. Solid squares are from 
direct measurements and blue triangles are the product of (Aa;|)/(Aa;j_) and t\\/t± (see text). The solid line is an empirical 
relation for visual guide. Inset illustrates schematically the distinct particle motions with positive (blue) and negative (red) 
translation-rotation couplings, (c) The extremum of C(t) v.s. D\\/D±. The solid line is a linear fit. The “prohibited zone” of 
Brownian motion is indicated by the red areas in (b) and (c). 



FIG. 4: (Color online). D\\/D± as a function of n for ellipsoids 
of different aspect ratios, p. Inset shows the linear relation 
between (1 — D\\/D ±)/(1 + D\\/D±_) and ( p 2 — 1)/( p 2 + 1), as 
predicted by the hydrodynamic model (Eq.[2J), where D\\/D_\_ 
are taken at n — 40no. 


same transition from the super-diffusive to the diffusive 
behaviors, allowing us to extract Di at long times. Con¬ 
sistent with the analysis of the lab-frame correlations, 
D\\/D±_ decreases with n and enters into “the prohibited 
zone” of Brownian diffusion at high n with D\\/D±_ < 1 
(Fig. (3 Jd) . The link between the coupling mode and the 
anisotropic body-frame diffusion is further illustrated by 
the linear relation between D\\/D± and the extremum 
of C(t ), where the emergence of negative C(t) coincides 
with D\\/D±_ < 1 (Fig. [3fc) . D\\/D±_ shows a qualita¬ 
tively similar trend for ellipsoids of other aspect ratios, 
p (Fig. |4|). But the variation of D\\/D± with n reduces 
with decreasing p. 

To understand the origin of the anomalous anisotropic 
diffusions, we apply the Green-Kubo formula (3l|, where 


Di is expressed as the integral of the velocity auto¬ 
correlation, Di = J 0 °° dt(vi(t + to)vi(to)), which is fur¬ 
ther approximated as Di ~ (v 2 ) J 0 °° dt(rii(t + to)rii(to)). 
Here, rii(t) = Vi(t)/\vi(t)\ is the direction of particle ve¬ 
locity vi along the particle’s major (i = ||) or minor 
(i = _L) axis at t. Thus, Di is decomposed into two 
parts: (1) the average particle velocity in the super- 
diffusive regime quantified by MSDs in a short-time inter¬ 
val At, (v 2 ) = (Ax 2 )/At 2 ; and (2) the auto-correlation 
of the velocity direction, Ai(t) = (rii(t + to)rii(to)). 
Aj( t) decays exponentially following Ai(t) = exp(— t/ri) 
[27[, where the correlation time, t*, indicates the per¬ 
sistence of the motion. Hence, Di = {vf)ri, leading to 
D\\/D± = ((A^)/(Ax 2 _))(t||/t j _) (Fig. Ed). At low n, 
while t || and r± are comparable (Fig. [SJd), (Ax|) is con¬ 
sistently larger than (Axj_) (Fig. [5^), a feature similar to 
Brownian motion of anisotropic particles [HI, [32]. How¬ 
ever, as n increases, (Ax^)/(Ax 2 ^} approaches to a con¬ 
stant close to 1, whereas r±_ grows larger than ry, result¬ 
ing in the anomalous D\\/D±_ < 1. Thus, the anomalous 
coupling arises from the more persistent motion of an 
anisotropic particle along its minor axis. 

To quantitatively interpret the anomalous diffusions 
from a first-principles calculation, one needs to con¬ 
sider hydrodynamic interactions among a large number 
of swarming bacteria, which is notoriously complicated 
and is still far from well understood M- Such a the¬ 
oretical investigation is out of the scope of our current 
experimental study. Nevertheless, we find that the de¬ 
creasing trend of T||/t_i_ emerges even when only consid¬ 
ering the far-field dipole flow of a single bacterium (see 
SI for detailed discussions Ezh 

In short, we extend previous models on enhanced diffu¬ 
sion of point-like spherical particles H 2 H 23 by assuming 
point-like particles with anisotropic shapes. The intrin- 
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FIG. 5: (Color online). Origin of the anomalous translation-rotation coupling, (a) Ratio of short-time MSDs along the major 
and minor axes in the super-diffusive regime with At = 0.1 s. (b) Ratio of correlation times along the major and minor axes. 
Schematics showing the average effective straining on ellipsoids induced by the extensile dipole flow of pushers such as E. coli 
(c) and the contractile dipole flow of pullers such as C. reinhardtii (d). 


sic correlation between the translation and rotation of an 
anisotropic point particle can be quantified using a scalar 
S = u • ( ), where v and co are the translational and 

\\axv\ 

angular velocity of the particle in a dipole flow created 
by a swimming bacterium [33|, and a is the unit vector 
along the major axis of the particle that forms an acute 
angle with v (Figs. [5b, d). We demonstrate that the spa¬ 
tial average of S,(S) , is negative for the extensile dipole 
flow of E. coli [27[. The negative sign indicates that the 
particle rotates in a direction that aligns its minor axis 
with the direction of its own translation, similar to the 
motion of a particle in a straining flow (Fig.[5h), whereas 
the magnitude, 1(5)1, gives the average speed of such a 
rotation. Under the influence of random fluctuations in 
a bacterial bath, this effective straining leads to (27| : 

T|| _ l + 2(S)t c /TT 

t_l i-2{S)t c y {> 

where the dimensionless quantity, 2| (S) |t c /7r, is the ratio 
of the correlation time of random fluctuations 0 , £ c , to 
the characteristic time of the straining-induced rotation, 
(tt/2)/\(S)\. Eq.H] successfully predicts < 1 for 

(S) < 0. Moreover, 2\{S)\t c /i: should be zero at n = 0 
and increase with n. Indeed, a linear approximation, 
2|(5)|£ c /7r = c(n/no), quantitatively fits our experiments 
with a numerical constant c = 1.6 x 10 3 (Fig. [SJd). Al¬ 
though the assumption of our current analysis only vali¬ 
dates at low n, where the far-field dipole flow dominates 
the bacteria-particle interaction, the quantitative agree¬ 
ment between experiments and this linear approximation 
indicates that at the coarse-grained level the swarming 
of a group of bacteria may still be treated as an effective 
dipole, albeit with a greatly enhanced dipolar strength. 
Last, the simple hydrodynamic model also predicts that 
the strength of the local straining decreases with the ellip¬ 
soids’ aspect ratio, p, and vanishes for spherical particles 


with p = 1, which, at high n, leads to a relation (27| : 

1 -(n/r±) 1 -(D\\/D±) ~ ( ) 

1 + (t||/tj_) ~ 1 + (D\\/D ± ) p 2 + 1’ 

consistent with our experiments (Fig. 0] inset). 

A self-propelled microswimmer can be generally cate¬ 
gorized as either a “pusher” or a “puller” depending on 
the far-field flow it creates HR Our simple calculation 
showed that the anomalous translation-rotation coupling 
is a consequence of extensile dipole flows—a defining fea¬ 
ture of pushers. Thus, this coupling should be universal 
for all pusher-type active fluids. In contrast, for puller 
microswimmers such as Chlamydomonas reinhardtii , (5) 
is positive (Fig. [5ji) [I?}. Therefore, an anisotropic par¬ 
ticle in puller-type of active fluids should exhibit an en¬ 
hanced positive translation-rotation coupling. However, 
due to the lack of the swarming behavior of pullers [2], the 
effect is likely to be weak. Our preliminary experiments 
with C. reinhardtii are consistent with this prediction. 

In addition to providing a new insight into the intrin¬ 
sic properties of active fluids, our results also deepen the 
understanding of transport processes in microbiological 
systems We have demonstrated that the transla¬ 

tion and rotation of anisotropic particles—naturally the 
most abundant particles in microbiological systems (e.g. 
dead bacterial bodies and macromolecules secreted by 
bacteria)—are not two decoupled degrees of freedom. In¬ 
stead, the rotation of an anisotropic particle profoundly 
influences how it translates and thus fundamentally mod¬ 
ifies its transport dynamics. 
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